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ABSTRACT 

We reconstruct the dark energy density px{z) as a free function from current type la supernova (SN 
la) data (Tonry et al. 2003; Barris et al. 2003; Knop et al. 2003), together with the Cosmic Microwave 
Background (CMB) shift parameter from CMB data (WMAP, CBI, and ACBAR), and the large scale 
structure (LSS) growth factor from 2dF galaxy survey data. We parametrize px{z) as a continuous 
function, given by interpolating its amplitudes at equally spaced z values in the rcdshift range covered 
by SN la data, and a constant at larger z (where px(z) is only weakly constrained by CMB data). We 
assume a flat universe, and use the Markov Chain Monte Carlo (MCMC) technique in our analysis. We 
find that the dark energy density px{z) is constant for < z < 0.5 and increases with redshift z for 
0-5 < z < 1 at 68.3% confidence level, but is consistent with a constant at 95% confidence level. For 
comparison, we also give constraints on a constant equation of state for the dark energy. 

Flux- averaging of SN la data is required to yield cosmological parameter constraints that are free of 
the bias induced by weak gravitational lensing (Wang 2000b). We set up a consistent framework for 
flux-averaging analysis of SN la data, based on Wang (2000b). We find that flux-averaging of SN la data 
leads to slightly lower f2 m and smaller time- variation in px(z). This suggests that a significant increase 
in the number of SNe la from deep SN surveys on a dedicated telescope (Wang 2000a) is needed to place 
a robust constraint on the time-dependence of the dark energy density. 

Subject headings: cosmology: observations - distance scale - supernovae:general 



1. INTRODUCTION 

Observational data of type la supernovae (SNe la) indi- 
cate that our universe is dominated by dark energy today 
(Riess et al. 1998; Perlmutter et al. 1999). The nature of 
dark energy is one of the great mysteries in cosmology at 
present. The time-dependence of the dark energy density 
px(z) can illuminate the nature of dark energy, and help 
differentiate among the various dark energy models (for ex- 
ample, Freese et al. (1987); Peebles & Ratra (1988); Frie- 
man et al. (1995); Caldwell, Dave, & Steinhardt (1998); 
Dodelson, Kaplinghat, & Stewart (2000); Deffayet (2001); 
Albrecht et al. (2002); Boyle, Caldwell, & Kamionkowski 
(2002); Freese & Lewis (2002); Griest (2002); Sahni & 
Shtanov (2002); Carroll, Hoffman, & Trodden (2003); Far- 
rar & Peebles (2003). See Padmanabhan (2003) and Pee- 
bles & Ratra (2003) for reviews with more complete lists 
of references) . 

SNe la can be calibrated to be good cosmological stan- 
dard candles, with small dispersions in their peak lumi- 
nosity (Phillips 1993; Riess, Press, & Kirshner 1995). The 
measurements of the distance-redshift relations of SNe la 
are most promising for constraining the time variation of 
the dark energy density px(z). The luminosity distance 
dh{z) = (X + z)r(z), with the comoving distance r(z) given 

by 

dz' 



state wx{z) = wq + w\z, then 



Px(z) 
Px(0) 



= 3u)iz _|_ ^3(l+ioo— iui) 
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e(z) = [n ro (i + z f + n fc (i + zf + n xPx (z)/ P 
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with fifc = 1 — £l m — Qx ■ If the dark energy equation of 
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The dark energy density p x (z) = Px(0)(l + z) 3 ( 1 + , "°) if 
the dark energy equation of state is a constant given by 
w . 

Most researchers have chosen to study dark energy by 
constraining the dark energy equation of state wx ■ How- 
ever, due to the smearing effect (Maor, Brustein, & Stein- 
hardt 2001) arising from the multiple integrals relating 
wx(z) to the luminosity distance of SNe la, it is 

extremely hard to constrain wx using SN data without 
making specific assumptions about wx (Barger & Marfa- 
tia 2001; Huterer & Turner 2001; Maor et al. 2002; Wasser- 
man 2002). If we constrain the dark energy density P x{z) 
instead, we minimize the smearing effect by removing one 
integral (Wang & Garnavich 2001; Tegmark 2002; Daly & 
Djorgovski 2003). 

It is important that there are a number of other probes 
of dark energy that arc complementary to SN la data (for 
example, see Podariu & Ratra (2001); Schulz & White 

(2001) ; Bean & Melchiorri (2002); Hu (2002); Sereno 

(2002) ; Bernstein & Jain (2003); Huterer & Ma (2003); 
Jimenez (2003); Majumdar & Mohr (2003); Mukherjee et 
al. (2003); Munshi & Wang (2003); Munshi, Porciani, & 
Wang (2003); Seo & Eisenstein (2003); Viel et al. (2003); 
Weller & Lewis (2003); Zhu & Fujimoto (2003)). Since 
different methods differ in systematic uncertainties, the 
comparison of them will allow for consistency checks, while 
the combination of them could yield tighter constraints on 
dark energy (for example, see Gerke & Efstathiou (2002); 
Hannestad & Mortsell (2002); Kujat et al. (2002)). 
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The most pressing question about dark energy that can 
be addressed by observational data is whether the dark 
energy density varies with time. In order to constrain the 
time- variation of dark energy density in a robust manner, 
it is important that we allow the dark energy density to 
be an arbitrary function of redshift z (Wang & Garnavich 
2001; Wang & Lovelace 2001; Wang et al. 2003). In this 
paper, we present a model-independent reconstruction of 
the dark energy density px{z), using SN la data published 
recently by the High-z Supernova Search Team (HZT) 
and the Supernova Cosmology Project (SCP) (Tonry et 
al. 2003; Barris et al. 2003; Knop et al. 2003), together 
with constraints from CMB (WMAP (Bennett et al. 2003), 
CBI(Pearson et al. 2003), and ACBAR(Kuo et al. 2002)) 
and large scale structure (LSS) data from the 2dF galaxy 
survey (Percival et al. 2002). 

Note that for clarity of presentation, we will label the 
samples of SNe la that we use according to the papers in 
which they were published. Hence we will refer to the 194 
SNe la from Tonry et al. (2003) and Barris et al. (2003) as 
the "Tonry/Barris sample" , and the 58 SNe la from Knop 
et al. (2003) as the "Knop sample" . 

Flux-averaging of SN la data is required to yield cosmo- 
logical parameter constraints that are free of the bias in- 
duced by weak gravitational lensing (Wang 2000b) . In this 
paper, we set up a consistent framework for flux-averaging 
analysis of SN la data, based on Wang (2000b). 

We assume a flat universe, and use the Markov Chain 
Monte Carlo (MCMC) technique in our analysis. 

Sec. 2 contains a consistent framework for flux-averaging 
analysis. We present our constraints on dark energy in 
Sec. 3. Sec. 4 contains a summary and discussions. 

2. A CONSISTENT FRAMEWORK FOR FLUX- AVERAGING 
ANALYSIS 

Since our universe is inhomogeneous in matter distribu- 
tion, weak gravitational lensing by galaxies is one of the 
main systematics 2 in the use of SNe la as cosmological 
standard candles (Kantowski, Vaughan, & Branch 1995; 
Frieman 1997; Wambsganss et al. 1997; Holz & Wald 1998; 
Metcalf & Silk 1999; Wang 1999; Wang, Holz, & Mun- 
shi 2002; Munshi & Wang 2003). Flux-averaging justifies 
the use of the distance-redshift relation for a smooth uni- 
verse in the analysis of type la supernova (SN la) data 
(Wang 2000b). Flux-averaging of SN la data is required 
to yield cosmological parameter constraints that are free 
of the bias induced by weak gravitational lensing (Wang 
2000b). 3 Here we set up a consistent framework for flux- 
averaging analysis of SN la data, based on Wang (2000b) . 



where F tr (z\s tr ) is the predicted flux due to the true cos- 
mological model parametrized by the set of cosmological 
parameters {s* r }, AF int is the uncertainty in SN la peak 
brightness due to intrinsic variations in SN la peak lumi- 
nosity and observational uncertainties, and /i is the magni- 
fication due to gravitational lensing by intervening matter. 
Therefore 
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The flux-averaging described in the Section 2.3 leads to 
the flux in each redshift bin 

F(z ibin ) = F tr {*iuJ (»)iun + {^F int ) iun . (7) 
For a sufficiently large number of SNe la in the z-th bin, 
(M)»nn=l- Hence 

N bin 

-E 



x Nbi M r ) 



N bin 



KnAF, 



int/ibi: 



[(AF, 



ibin 



'F,i bi . 



'F,i bi . 



(8) 

Comparison of Eq.(8) and Eq.(6) shows that flux- 
averaging can remove/reduce the gravitational lensing ef- 
fect, and leads to a smaller y 2 per degree of freedom for the 
true model, compared to that from without flux-averaging. 

2.2. Flux statistics versus magnitude statistics 

Normally distributed measurement errors are required if 
the x 2 parameter estimate is to be a maximum likelihood 
estimator (Press et al. 1994). Hence, it is important that 
we use the \ 2 statistics with an observable that has a error 
distribution closest to Gaussian. 

So far, it has been assumed that the distribution of ob- 
served SN la peak brightness is Gaussian in magnitudes. 
Therefore, for a given set of cosmological parameters {s} 

,2 _ Mgj) - t4)( z *\ S )f 



X 



E 



(9) 



2.1. Why flux- averaging ? 

The reason that flux-averaging can remove/reduce grav- 
itational lensing bias is that due to flux conservation, the 
average magnification of a sufficient number of standard 
candles at the same redshift is one. 

The observed flux from a SN la can be written as 



F(z) = F mt n, 



Fint = F tr (z\s tr ) + AFi 



in!, • 



(4) 



where ^{z) — 5 log (c?L(z)/Mpc) + 25, and c?l(z) = 
(1 + z)r(z) is the luminosity distance. 

However, while we do not have a very clear understand- 
ing of how the intrinsic dispersions in SN la peak luminos- 
ity is distributed, the distribution of observational uncer- 
tainties in SN la peak brightness is Gaussian in flux, since 
CCD's have replaced photometric plates as detectors of 
photons. 



2 The others systematics are possible gray dust (Aguirre 1999) and SN la peak luminosity evolution (Drell, Loredo, & Wasserman 2000; Riess 
et al. 1999; Wang 2000b); so far, there is no clear evidence of either. 

3 To avoid missing the faint end (which is fortunately steep) of the magnification distribution of observed SNe la, only SNe la detected well 



above the threshold should be used in flux-averaging. 
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In this paper, we assume that the intrinsic dispersions 
in SN la peak brightness is Gaussian in flux, and not in 
magnitude as assumed in all previous publications. The 
justifications for this preference will be presented in detail 
elsewhere (Wang et al. 2004) 4 . Thus, 



[F{ Zi )-Fv{ Zi \s)Y 



(10) 



'FA 



Since the peak brightness of SNe la have been given in 
magnitudes with symmetric error bars, m pea k ± a m , we 
obtain equivalent errors in flux as follows: 

_ F(m peak + er m ) - F(m peak - a. m ) 
g f = ^ ^ U > 

We will refer to Eq.(9) as "magnitude statistics", and 
Eq.(10) as "flux statistics". For reference and comparison, 
we present results in both "magnitude statistics" and "flux 
statistics" in this paper. However, a consistent framework 
for flux-averaging is only straightforward in "flux statis- 
tics". 5 

2.3. A recipe for flux-averaging 

The procedure for flux-averaging in Wang (2000b) is 
for minimizing y 2 using the subroutines from Numeri- 
cal Recipes (Press et al. 1994). As described in Wang 
(2000b), the fluxes of SNe la in a redshift bin should 
only be averaged after removing their redshift dependence, 
which is a model-dependent process. For x 2 statistics us- 
ing MCMC or a grid of parameters, here are the steps in 
flux-averaging: 

(1) Convert the distance modulus of SNe la into 
"fluxes", 

-2 



F(Zj) = IQ-^o^i)- 25 )/ 2 - 5 = 



jdata 



Mpc 



(12) 



(13) 



(2) For a given set of cosmological parameters {s}, ob- 
tain "absolute luminosities", {C(zj)}, by removing the 
redshift dependence of the "fluxes", i.e., 

£(z,) = d|(z»F(z 3 ). 

(3) Flux-average the "absolute luminosities" 

each redshift bin i to obtain j : 

-i 1 ^ , 1 



N 



3=1 



N 

H z r ( 14 ) 
of the z-th redshift 



(4) Place C at the mean redshift Zt 
bin, now the binned flux is 

F(z i )=r/dl(z i \s). (15) 

The l-o- error on each binned data point F l , af , is taken to 
be the root mean square of the l-o errors on the unbinncd 
data points in the i-th redshift bin, {F*} (j = 1, 2, N), 

multiplied by 1/VN (see Wang (2000a)). 

(5) For the flux-averaged data, {F(z~i)}, we find 

[F(zj) -FP(^|s)]" 



X 



= E 



(16) 



3. CONSTRAINTS ON DARK ENERGY 

The Tonry/Barris SN la sample consists of 194 SNe la 
with z > 0.01 and extinction A v < 0.5 (Tonry et al. 2003; 
Barris et al. 2003). To examine the effect of the two SNe 
la at the high reshift end (z — 1.199 and z = 1.755 respec- 
tively), we also present the results for 193 SNe la (omitting 
the SN la at z = 1.755) and 192 SNe la (omitting the two 
SNe la at z = 1.199 and z = 1.755). 

The Knop SN la sample consists of 58 SNe la (the "All 
SCP SNe" data set from Knop et al. (2003)). These data 

should be compared with m e J^ = 51og(7J ^i) + offset, 
with 



offset = 5 log (2997.9//i) + 25 + M s 



(17) 



where M$n is the peak absolute magnitude of SNe la. 
Note that for the Knop sample, the flux statistics must be 
done with a revised definition of flux, 



F(z) 



[H d L (z)} 
Mpc 



data 



10 



-2(mf f/ -0ffset)/5 



(18) 

to be compared with theoretical prediction of F p (z\s) = 
[Hq di(2;|s)/Mpc]~ 2 for a given set of cosmological param- 
eters {s}. 

Note that an additional uncertainty from the redshift 
dispersion due to peculiar velocity must be added to the 
uncertainty of each SN la data point. The Knop sam- 
ple already include a dispersion of 300 km/s along the 
line of site. To add 500 km/s dispersion in z to the 
SN data in the Tonry/Barris sample, one must propagate 
o z = c _1 500 km/s into an additional uncertainty in the 
luminosity distance cIl{z), then add it in quadrature to 
published uncertainty in oIl(z). Note that this process is 
dependent on the cosmological model, and must be done 
for each set of cosmological parameters during the likeli- 
hood analysis (Riess et al. 1998; Wang 2000b). 

To obtain tighter constraints on dark energy, we also in- 
clude constraints from CMB (WMAP, CBI, ACBAR) and 
LSS (2dF) in our analysis. Since CMB data clearly indi- 
cate that we live in a flat universe, we assume D, m +Qx = 1 
in all our results. 6 

We use the Markov Chain Monte Carlo (MCMC) tech- 
nique (see Neil (1993) for a review), illustrated for exam- 
ple in Lewis & Bridle (2002), in the likelihood analysis. 
At its best, the MCMC method scales approximately lin- 
early in computation time with the number of parameters. 
The method samples from the full posterior distribution 
of the parameters, and from these samples the marginal- 
ized posterior distributions of the parameters can be es- 
timated. We have derived all our probability distribution 
functions (pdf) of the cosmological parameters from 10 6 
MCMC samples. 



where (^|si) = (d L (z\s) /Mpc) . 

4 Our study of intrinsic peak luminosities of nearby SNe la shows that their distribution is much more Gaussian in flux than in magnitude. 

5 If the the dispersions in SN la peak brightness were Gaussian in magnitude, flux-averaging would introduce a small bias. 

6 Allowing both S7 m and !2x to vary would lead to greatly increased uncertainty in dark energy constraints, such that no interesting constraints 
can be obtained from current data. 
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3.1. The likelihood analysis 



We use a % 2 statistic 



2 _ 2 | 2 | 2 

X — XSN + XCMB + XiSSi 



(19) 

where Xsw * s gi ven by Eq.(16) and (10) for flux statis- 
tics (with and without flux-averaging), and Eq.(9) for 
magnitude statistics. Xcmb an d Xlss are contributions 
from CMB and LSS data respectively. The likelihood 
C oc e~ x I 2 if the measurement errors are Gaussian (Press 
et al. 1994). 

When the cosmological parameters are varied, the shift 
in the whole CMB angular spectrum is determined by the 
shift parameter (Bond, Efstathiou, & Tcgmark 1997; Mcl- 
chiorri et al. 2002; Odman at el. 2002) 



U = H r(z dec ) (20) 

where r(zdec) denotes the comoving distance to the decou- 
pling surface in a flat universe. Note that this is a robust 
way to include CMB constraints since the CMB depends 
on f2 m and h in the combination of the physical parame- 
ter n m h 2 . The results from CMB (WMAP, CBI, ACBAR) 
data correspond to 1Z = 1.716 ± 0.062 (using results in 
Spergel et al. (2003)). We include the CMB data in our 

analysis by adding Xcmb = [(^ — ^o) / ' > where 1Z is 
computed for each model using Eq.(20). 

Following Knop et al. (2003), we include the LSS con- 
straints from 2dF in terms of the growth parameter / = 
dhiD I cflna, where a is the cosmic scale factor, and D is the 
linear fluctuation growth factor, D{t) = 5^(x,t)/5(x), 
given by 

D(t) + 2H(z)D(t) - | il m H 2 (l + zfD{t) = 0, (21) 

where the dots denote derivatives with respect to t. 

The Hubble parameter H{z) = H E{z) (see Eq.(2)). 
Since (3 — f/bi, the 2dF constraints of (3(z ~ 0.15) = 
0.49 ± 0.09 (Hawkins et al. 2003) and 6i = 1.04 ± 0.11 
(Verde et al. 2002) yields / = f (z = 0.15) = 0.51 ± 0.11. 
We include the 2dF constraints in our analysis by adding 

Xlss = {lf( z ~ 0-15) — fo]/o-f } , where / = dhiD/dhia 
is computed for each model using D obtained by numeri- 
cally integrating Eq.(21). 

Note that we have chosen to use only the most conserva- 
tive and robust information, the CMB shift parameter and 
the LSS growth factor, from CMB and LSS observations 7 . 
It is important that the limits on these are independent 
of the assumption on dark energy made in the CMB and 
LSS data analysis. Further, by limiting the amount of in- 
formation that we use from CMB and LSS observations 
to complement the SN la data, we minimize the effect of 
the systematics inherent in the CMB and LSS data on our 
results. 

3.2. Constraints on a constant dark energy equation of 

state wo 

The most popular and simplest assumption about dark 
energy is that it has a constant equation of state wq . Here 
we present constraints on a constant dark energy equation 
of state. 

Fig. 1 shows the marginalized pdf of the matter density 
fraction Cl m , the dimensionless Hubble constant h, and the 
constant dark energy equation of state wq . 

7 These observations provide a vast amount of information as detailed 



The first four rows of figures in Fig.l are results ob- 
tained using the Tonry/Barris SN la sample, requiring that 
z > 0.01 and extinction Ay < 0.5 (which yields a total of 
194 SNe la). To examine the effect of two SNe la at the 
high reshift end (z = 1.199 and z — 1.755 respectively), 
we also present the results for 193 SNe la (omitting the 
SN la at z = 1.755) and 192 SNe la (omitting the two SNe 
la at z = 1.199 and z — 1.755). The solid, dotted, and 
dashed lines indicate the results for 192, 193, and 194 SNe 
la respectively. 

The first two rows of figures in Fig.l are results for SN 
la data from the Tonry/Barris sample only, without (first 
row) and with (second row) flux-averaging (Az — 0.05). 
Note that inclusion of the two highest redshift SNe la at 
z = 1.199 and z = 1.755 leads slightly higher O m and more 
negative w - Flux-averaging leads to broader pdf for Q m , 
with somewhat lower mean TO , and broader pdf for h. 

The 3rd and 4th rows of figures in Fig.l are results for 
SN la data from the Tonry/Barris sample, combined with 
constraints from CMB (WMAP, CBI, ACBAR) and LSS 
(2dF) data. The inclusion of the two highest redshift SNe 
la at z — 1.199 and z = 1.755 makes less difference in 
the estimated parameters, since the inclusion of the CMB 
and LSS data reduces the relative weight of these two data 
points. The main effect of flux-averaging is a broader pdf 
for h. 

The 5th row of figures in Fig.l are results obtained us- 
ing 58 SNe la from the Knop sample, using flux-averaged 
statistics (solid) and magnitude statistics (dotted) respec- 
tively. Note flux averaging significantly broadens all the 
pdf's. The central figure is equivalent to a pdf in h (see 
Eq.(17). These are consistent with similar results derived 
using the SNe la from the Tonry/Barris sample (3rd row 
of figures in Fig.l). 

Table 1 gives the marginalized 68.3% and 95% confi- 
dence level (C.L.) of fl m , h, and wq. These have been 
computed using 10 6 MCMC samples. 



3.3. Constraints on dark energy density as a free 
function 

To place model-independent constraints on dark energy, 
we parametrize px(z) as a continuous function, given by 
interpolating its amplitudes at equally spaced z values in 
the redshift range covered by SN la data (0 < z < z max ), 
and a constant at larger z (z > z max , where px{z) is 
only weakly constrained by CMB data). The values of 
the dimensionless dark energy density /j = p(zi)/px(0) 
(i = 1,2, ...,rif) are the independent variables to be esti- 
mated from data. We interpolate px{z) using a polyno- 
mial of order n/ for < z < z max . 

Since the present data can not constrain px(z) for 
nj > 2, we present results for rif — 2, i.e., with px{z) 
parametrized by its values at z = z maz /2, z maz - 

Fig. 2 shows the marginalized pdf of the matter density 
fraction TO (column 1), the dimensionless Hubble con- 
stant h (column 2), and dimensionless dark energy density 
at z = z max /2 and z = z max (columns 3 and 4), obtained 
using current SN la data (Tonry et al. 2003; Barris et al. 
2003; Knop et al. 2003), flux-averaged and combined with 



in the publications from the WMAP and 2dF teams. 
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CMB (WMAP, CBI, ACBAR) and 2dF data. The first 
three rows of figures are results for 192, 193, and 194 SNe 
la from the Tonry/Barris sample, while the fourth row are 
results for 58 SNe la from the Knop sample. The dark 
energy density at z = z max is not well constrained when 
the z = 1.755 SN la is included in the analysis; this is as 
expected since this extends p x (z) to z max — 1.755, with 
only one SN la at z > 1.2. 

Note that when the estimated parameters are well con- 
strained, flux averaging generally leads to slightly lower 
estimates of VL m and p X (z) (at z — z max /2 and z — z max ). 

Table 2 gives the marginalized 68.3% and 95% C.L. of 
fi m , h, px(z max /2)/p x (0), and px{z ma x)l Px(0)- These 
have been computed using 10 6 MCMC samples. 

Fig. 3 shows the dark energy density px{z) recon- 
structed from current SN la (the Tonry/Barris sample 
and the Knop sample), CMB (WMAP, CBI, ACBAR) 
and LSS (2dF) data. The heavy (light) lines indicate 
the 68.3% (95%) C.L. of the reconstructed px(z). The 
dot-dashed line indicates the cosmological constant model, 
px{z)/px(0) = 1. The 68.3% and 95% C.L.'s of p x (z) 
are marginalized confidence levels, computed at each z 
using 10 6 MCMC samples, with the correlation between 
Px{-5z m ax) and p x {z ma x) fully included. 

Fig. 3(a) shows the reconstructed px(z) using 192, 
193, and 194 SNe la from the Tonry/Barris sample, 
flux-averaged and combined with CMB (WMAP, CBI, 
ACBAR) and LSS (2dF) data. The densely (sparsely) 
shaded regions are the 68.3% (95%) C.L. of p x (z) for 192 
SNe la (at z < 1.056). The heavy (light) dotted and 
dashed lines are the 68.3% (95%) C.L. of p x (z) for 193 
and 194 SNe la respectively. Note that the px(z) recon- 
structed from 193 SNe la (adding the SN la at z = 1.199) 
nearly overlaps from that from 192 SNe la for z < 1.056. 
However, the px{z) reconstructed from 194 SNe la (adding 
the SNe la at z = 1.199 and z = 1.755) deviates no- 
tably from that from 192 SNe la for 0.7 < z < 1.056, 
although the 68.3% C.L. regions overlap. Clearly, the re- 
constructed px{z) is constant for < z < 0.5 and in- 
creases with redshift z for 0.5 < z < 1 at 68.3% C.L., 
but is consistent with a constant at 95% C.L. We note 
that at 90% C.L., p x {.5 z max )/p(0) = [0.83,1.59], and 
Px(zmax)/p(0) = [1.03,6.85]; this indicates that px(z) 
varies with time at approximately 90% C.L. 

Fig. 3(b) shows the reconstructed px{z) using 192 SNe 
la from the Tonry/Barris sample (same as in Fig. 3(a)) and 
that from the 58 SNe la of the Knop sample (dotted lines). 
The 68.3% C.L. regions overlap. However, the 58 SCP SNe 
la seem to favor px{z) < 1 at 0.5 < z < 1, and has much 
larger uncertainties at z > 0.5. 

3.4. Comparison with previous work 

The SN observational teams have published their data 
together with constraints on a constant dark energy equa- 
tion of state (Tonry et al. 2003; Knop ct al. 2003). These 
results should be compared with our results for a constant 
wx{z) using magnitude statistics (see Table 1). Using a 
2dF prior of Q m h = 0.20 ± 0.03 (Percival et al. 2002) 
and assuming a flat universe, Tonry et al. (2003) found 
that -1.48 < w < -0.72 at 95% C.L.; this is close to 
-1.24 < w < -0.74 at 95% C.L. (using 192 SNe la 
together with CMB and LSS constraints as discussed in 



Sec.3.1) from Table 1. Using the same CMB and LSS 
constraints as us and assuming a flat universe, Knop et 
al. (2003) found that w = -1.05±g;^ (statistical) ±0.09 
(identified systematics) at 68.3% C.L.; close to our results 
of w = -0.99 ± 0.16 at 68.3% C.L. 

At the completion of our analysis, we became aware that 
Alam ct al. (2003) and Choudhury & Padmanabhan (2003) 
have found that current SN la data favor wx (z) < —1. Al- 
though our results are qualitatively consistent with these, 
there are significant differences in both analysis technique 
and quantitative results. 

Alam ct al. (2003) used 172 SNe la from Tonry et al. 
(2003), and obtained a reconstructed wx(z) that deviates 
quite significantly from wx{z) = — 1. In their paper (v2), 
Figs. 3 and 14 (both assuming TO = 0.3) arc consistent 
with what we found. We note that some of their recon- 
structed Wx{z)'s (see their Figs. 4, 6, 8, 10, and 16) have 
decreasing errors for z > 1.2 (where there are only two SNe 
la). This illustrates the fact that while free fitting forms 
for px(z) or Wx{z) generally give increasing errors with 
redshift (large error where there are very few observed SNe 
la, see Fig. 3 of this paper), this may not be true for other 
parametrizations (such as used in Alam et al. (2003)). 

Choudhury & Padmanabhan (2003) used 194 SNe la 
(the Tonry/Barris sample), and presented their results for 
w x (z) = w + W! z/(l + z) with n m = 0.29, 0.34, and 0.39 
(no marginalization over TO ). 

4. SUMMARY AND DISCUSSION 

In order to place model-independent constraints on dark 
energy, we have reconstructed the dark energy density 
px(z) as a free function from current SN la data (Tonry 
et al. 2003; Barris et al. 2003; Knop ct al. 2003), to- 
gether with CMB (WMAP, CBI, and ACBAR) and LSS 
(2dF) data. We find that the dark energy density px(z) 
is constant for < z < 0.5 and increases with redshift z 
for 0.5 < z < 1 at 68.3% C.L., but is consistent with a 
constant at 95% C.L. (see Fig. 3). 

Flux-averaging of SN la data is required to yield cos- 
mological parameter constraints that arc free of the bias 
induced by weak gravitational lensing (Wang 2000b). We 
have developed a consistent framework for flux-averaging 
analysis of SN la data, and applied it to current SN la 
data. We find that flux-averaging of SN la data generally 
leads to slightly lower Cl m and smaller time-variation in 
Px(z). 

We note that flux-averaging of SNe la has more effect on 
the Knop sample than the Tonry/Barris sample. This may 
be due to the fact that the measurement errors of the ma- 
jority of the SNe la in the Tonry/Barris sample have been 
"Gaussianized" in magnitudes by averaging over several 
different analysis techniques (Tonry et al. 2003). How- 
ever, it is likely that SN la peak brightness distribution 
is Gaussian in flux, instead of magnitudes (Wang et al. 
2004). A consistent framework for flux-averaging is only 
straightforward (as presented in Sec. 2) if the distribution 
of SN la peak brightnesses is Gaussian in flux. Our re- 
sults suggest that observers should publish observed SN 
la peak brightnesses with uncertainties in flux, to allow 
detailed flux-averaging studies. 

Our results include an estimate of the Hubble constant 
H = h 100km/sMpc _1 from the Tonry/Barris sample of 
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194 SNe la. Since the Tonry/Barris sample data used a 
fixed value of hfi x — 0.65 (Tonry et al. 2003) in the de- 
rived distances, we divide their derived distances H (Il{z) 
by hfi x and marginalize over Hq in our analysis. Our 
MCMC method yields smooth pdf's for all marginalized 
parameters. The errors on the estimated h in Tables 1 & 
2 are statistical errors only, not including a much larger 
systematic error contributed by the intrinsic dispersion in 
SN la peak luminosity of <r™* ~ 0.17 magnitudes Hamuy 
et al. (1996). This implies a systematic uncertainty in h of 
7.83%, or erf* ~ 0.05 for the h values tabulated in Table 1 
& 2. This yields an estimate for h that overlap with those 
from Branch (1998) and Freedman et al. (2001) within 1 a. 

It is intriguing that the current SN la data, together 
with CMB and galaxy survey data, indicate that px{z) 



varies with time at approximately 90% C.L. (see Fig. 3 
and Sec. 3. 3). If the trend in px(z) that we have found 
is confirmed by future observational data, it will have rev- 
olutionary implications for particle physics and cosmology. 
Since the uncertainty in px{z) is large where there are few 
observed SNe la (see Fig. 3), we expect that a significant 
increase in the number of SNe la, obtained from dedi- 
cated deep SN la searches (Wang 2000a), will allow us to 
place robust and more stringent constraints on the time- 
dependence of the dark energy density. 

This work is supported in part by NSF CAREER grant 
AST-0094335. We are grateful to David Branch and 
Michael Vogeley for helpful comments. 



REFERENCES 



Aguirrc, A.N. 1999, ApJ, 512, L19 

Alam, U.; Salmi, V.; Saini, T.D.; & Starobinsky, A. A. 2003, astro- 
ph/0311364 

Albrccht, A.; Burgess, CP.; Ravndal, F.; & Skordis, C. 2002, Phys. 

Rev. D65, 123507 
Armendariz-Picon, C, Mukhanov, V., Steinhardt, P.J. 2000, Phys. 

Rev. Lett. 85, 4438 
Barger, V., & Marfatia, D. 2001, Phys. Lett. B, 498, 67 
Barris, B.J., et al. 2003, astro-ph/0310843, ApJ, in press 
Bean, R.; & Mclchiorri, A. 2002, Phys. Rev. D65, 041302 
Bennett, C, et al. 2003, ApJ, Suppl. 148, 1 

Bernstein, G.M.; & Jain, B. 2003, astro-ph/0309332, ApJ, in press 
Bond.J.R.; Efstathiou, G.; & Tegmark, M. 1997, MNRAS, 291, L33 
Boyle, L.A.; Caldwell, R.R.; Kamionkowski, M. 2002, Phys. Lett. 

B545, 17 
Branch, D. 1998, ARA&A, 36, 17 

Caldwell, R., Dave, R., Steinhardt, P. 1998, Phys. Rev. Lett., 80, 
1582 

Carroll, S.M.; Hoffman, M., & Trodden, M. 2003, Phys. Rev. D68, 
023509 

Choudhury, T.R., & Padmanabhan, T. 2003, astro-ph/0311622 
Daly, R.A., & Djorgovski, S.G. 2003, ApJ, 597, 9 
Deffayet, C. 2001, Phys. Lett. B502, 199 

Dodelson, S., Kaplinghat, M., & Stewart, E. 2000, Phys. Rev. Lett. 
85, 5276 

Drell, P.S.; Loredo, T.J.; & Wasserman, I. 2000, ApJ, 530, 593 
Farrar, G.R., & Peebles, P. J. E. 2003, astro-ph/0307316 
Freedman, W. L. et al. 2001, ApJ, 553, 47 

Freese, K., Adams, F.C., Frieman, J. A., and Mottola, E. 1987, Nucl. 

Phys. B287, 797 
Freese, K., and Lewis, M., 2002, Phys. Lett., B540, 1 
Frieman, J., Hill, J., Stebbins, A., and Waga, I. 1995, Phys. Rev. 

Lett., 75, 2077 
Frieman, J. A. 1997, Comments Astrophys., 18, 323 
Gcrke, B.F., Efstathiou, G. 2002, MNRAS, 335, 33 
Griest, K. 2002, Phys. Rev. D66, 123501 
Hamuy, M., et al. 1996, AJ, 112, 2408 

Hannestad, S., & Mortsell, E. 2002, Phys. Rev. D66, 063508 
Hawkins, E. et al. 2003, astro-ph/0212375, MNRAS in press 
Holz, D.E. & Wald, R.M. 1998, Phys. Rev., D58, 063501 
Hu, W. 2002, Phys.Rev. D66, 08351 
Huterer, D., Ma , C. 2003, astro-ph/0307301 
Huterer, D., & Turner, M.S. Phys. Rev. D64, 123527 
Jimenez, R. 2003, New Astron. Rev. 47, 761 
Kantowski, R.; Vaughan, T.; Branch, D. 1995, ApJ, 447, 35 
Knop, R. A., ct al. 2003, astro-ph/0309368, ApJ, in press 
Kujat, J.; Linn, A.M.; Scherrer, R.J.; & Weinberg, D.H. 2002, ApJ, 
572, 1 

Kuo, C.L., et al. 2002, submitted to ApJ, astro-ph/0212289 
Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511, astro- 
ph/0205436 

Maor, I., Brustein, R., & Steinhardt, P.J. 2001, Phys. Rev. Lett., 86, 
6; Erratum-ibid. 87 (2001) 049901 



Maor, I.; Brustein, R.; McMahon, J.; & Steinhardt, P.J. 2002, Phys. 

Rev. D65, 123003 
Mclchiorri, A.; Mersini, L.; Odman, C.J.; & Trodden, M. 2002, astro- 

ph/0211522 
Mctcalf, R. B., & Silk, J. 1999, ApJ, 519, LI 

Majumdar, S., & Mohr, J.J. 2003, astro-ph/0305341, submitted to 
ApJ 

Mukherjee, P.; Banday, A. J.; Riazuelo, A.; Gorski, K.M.; & Ratra, 

B. 2003, astro-ph/0306147, ApJ, in press 
Munshi, D.; Porciani, C; & Wang, Y. 2003, astro-ph/0302510, 

MNRAS, in press 
Munshi, D., and Wang, Y. 2003, ApJ, 583, 566 

Neil, R.M. 1993, ftp://ftp.cs.utoronto.ca/pub/ radford/rcview.ps.gz 
Odman, C.J.; Melchiorri, A.; Hobson, M.P.; & Lasenby, A.N. 2002, 

astro-ph/0207286 
Padmanabhan, T. 2003, Physics Reports 380, 235-320 
Peebles, P. J. E.; Ratra, B. 1988, ApJ, 325L, 17 
Peebles, P. J. E.; Ratra, B. 2003, Rev.Mod.Phys. 75, 559-606 
Pearson, T.J., et al. 2003, ApJ, 591, 556-574 
Percival, W.J. et al. 2002, MNRAS, 337, 1068 
Perlmutter, S., et al. 1999, ApJ, 517, 565 
Phillips, M.M., ApJ, 413, L105 (1993) 
Podariu, S.; & Ratra, B. 2001, ApJ, 563, 28 

Press, W.H., Teukolsky, S.A., Vettering, W.T., & Flannery, B.P. 

1994, Numerical Recipes, Cambridge University Press, Cambridge. 
Riess, A. G., ct al. 1998, AJ, 116, 1009 

Riess, A.G., Press, W.H., and Kirshner, R.P., ApJ, 438, L17 (1995) 

Riess, A.G., et al. 1999, AJ, 118, 2675 

Sahni, V., & Shtanov, Y. 2002, astro-ph/0202346 

Schulz, A.E.; & White, M. 2001, Phys. Rev. D64, 043514 

Seo, H.; & Eisenstein, D.J. 2003, astro-ph/0307460 

Sereno, M. 2002, Astron. Astrophys., 393, 757 

Spergel, D.N., et al. 2003, astro-ph/0302209 

Tegmark, M. 2002, Phys. Rev. D66, 103507 

Tonry, J.L., et al. 2003, ApJ, 594, 1-24 

Verde, L. et al. 2002, MNRAS, 335, 432 

Viel, M.; Matarrese, S.; Theuns, T.; Munshi, D.; & Wang, Y. 2003, 

MNRAS, 340, L47 
Wambsganss, J., Cen, R., Xu, G., & Ostriker, J. P. 1997, ApJ, 475, 

L81 

Wang, Y. 1999, ApJ, 525, 651 

Wang, Y. 2000a, ApJ, 531, 676 

Wang, Y. 2000b, ApJ, 536, 531 

Wang, Y., and Garnavich, P. 2001, ApJ, 552, 445 

Wang, Y., and Lovelace, G. 2001, ApJ, 562, L115 

Wang, Y., Holz, D.E., and Munshi, D. 2002, ApJ, 572, L15 

Wang, Y.; Freese, K.; Gondolo, P.; & Lewis, M. 2003, ApJ, 594, 25 

Wang, Y., ct al. 2004, in preparation 

Wasserman, I. 2002, Phys. Rev. D66, 123511 

Weller, J., & Lewis, A.M. 2003, MNRAS, in press 

Zhu, Z.; & Fujimoto, M. 2003, astro-ph/0312022, ApJ in press 



Table 1 

Estimated cosmological parameters (mean, 68.3% C.L., 95% C.L.) assuming wx(z) = wq 
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Tonry/Barris sample SNe 



192 SNe (z mo *=l-056) 

flux, binned 6 .47 
flux, unbinned .47 
mag.; unbinned .49 


.36, .57 
.43, .52 
.41, .57 


.14, .63 
.34, .55 
.22, .61 


.661 
.656 
.662 


h a 

.645, .674 
.644, .666 
.651, .671 


.631, .690 
.636, .676 
.644, .681 


w X 2 miJ N dof 

-2.37 [-3.56, -1.20] [-5.34, -.69] 13.34/19 
-3.08 [-4.09, -2.08] [-5.50, -1.46] 209.10/189 
-2.25 [-3.18, -1.34] [-4.36, -.80] 193.36/189 


193 SNe (z max =1.199) 

flux, binned 6 .48 
flux, unbinned .48 
mag.; unbinned .51 


.39, .58 
.43, .52 
.44, .58 


.18, .63 
.35, .56 
.29, .62 


.661 
.656 
.663 


.645, .676 
.645, .665 
.653, .671 


.631, .690 
.637, .676 
.644, .683 


-2.53 [-3.85, -1.30] [-5.54, -.73] 14.36/20 
-3.13 [-4.10, -2.13] [-5.64, -1.50] 210.56/190 
-2.47 [-3.42, -1.49] [-5.14, -.91] 194.86/190 


194 SNe (z moa: =1.755) 

flux, binned 6 .49 [.40, .58] [.20, .63] .661 [.645, .675] [.632, .690] -2.54 [-3.80, -1.34] [-5.62, -.76] 14.40/21 
flux, unbinned .48 [.44, .53] [.37, .56] .656 [.645, .665] [.637, .675] -3.17 [-4.13, -2.22] [-5.47, -1.58] 210.74/191 
mag.; unbinned .51 [.45, .58] [.30, .62] .663 [.653, .671] [.644, .683] -2.51 [-3.46, -1.55] [-4.97, -.95] 194.88/191 


192 SNe (z moa: =1.056) + CMB & LSS 
flux, binned 6 .28 
flux, unbinned .26 
mag.; unbinned .29 


.23, .33 
.22, .31 
.24, .34 


.19, .39 
.18, .36 
.21, .39 


.652 
.643 
.654 


.638, .665 
.634, .650 
.646, .661 


.627, .677 
.627, .657 
.639, .668 


-.95 [-1.09, -.82] [-1.27, -.72] 15.44/21 
-1.15 [-1.29, -1.00][-1.53, -.90] 216.76/191 
-.95 [-1.07, -0.83] -1.24, -.74] 196.60/191 


193 SNe (z mttx = 1.199) + CMB & LSS 

flux, binned 6 .29 [.24, .34] [.20, .39] .652 
flux, unbinned .27 [.22, 31] [.18, .37] .643 
mag.; unbinned .30 [.25, .35] [.21, .40] .654 


.638, .663 
.635, .650 
.645, .660 


.626, .676] -.95 
.627, .658] -1.15 
.638, .667] -.95 


-1.08, -.82] [-1.27, -.71] 17.04/22 
-1.29, -1.00] [-1.54, -.90] 219.00/192 
-1.07, -.82] [-1.25, -.74] 199.04/192 


194 SNe (z maa; =1.755) + CMB & LSS 

flux, binned 6 .29 [.24, .34] [.20, .40] .651 [.638, .663] [.625, .676] -.95 [-1.08, -.81] [-1.28, -.71] 17.44/23 
flux, unbinned .27 [.22, .32] [.18, .37] .642 [.634, .650] [.627, .657] -1.15 [-1.31, -1.00] [-1.57, -.90] 219.8/193 
mag.; unbinned .30 [.25, .35] [.21, .41] .654 [.645, .660] [.638, .668] -.95 [-1.07, -.82] [-1.27, -.73] 199.50/193 


58 Knop sample SNe (z max =0.863) 
+ CMB & LSS Q m offset-23.5 w *LiJ N lof 


flux, binned 6 .22 
flux, unbinned .21 
mag.; unbinned .27 


.15, .28 
.16, .26 
.21, .32 


.11, .36 
.13, .32 
.17, .39 


.414 
.395 
.367 


.347, .480 
.356, .434 
.331, .406 


.286, .547 
.321, .471 
.295, .441 


-1.20 
-1.18 
-.99 


[-1.45, -.96] [-1.73, -.74] 9.1/13 
-1.34, -1.02] [-1.54, -.88] 61.34/57 
-1.15, -.83] [-1.34, -.70] 55.86/57 



a Statistical error only, not including the contribution from the much larger SN la absolute magnitude error of a™ 1 ~ 0.05 (see Sec. 4). 
6 Flux-averaged with Az = .05. c The number of degrees of freedom. 

Table 2 

Estimated cosmological parameters (mean, 68.3% C.L., 95% C.L.) for arbitrary px(z), from SN la data combined with CMB and LSS data 











h a 




px (-5 z max ) / p x (0) 


Px{Zmax)/pX (0) 


~&min 1 ^ do f 


192 SNe 
binned flux 6 
unbinned flux 
unbinned mag. 


.33 
.34 
.34 


(Zmax =1.056) 

[.27,. 39] [.22, .46] 
[.28, .39] [.24, .45] 
[.28, .40] [.23, .46] 


Tonry/Barris sample 

.660 [.644, .673] [.630, .688] 
.655 [.645, .663] [.637, .671] 
.662 [.652, .670] [.645, .678] 


1.19 [.97, 1.42] [.76, 1.67] 
1.09 [.88, 1.31] [.67, 1.56] 
1.22 [1.01, 1.43][.80, 1.66] 


3.61 [1.84, 5.41] [.73, 7.53] 
5.02 [3.27, 6.82][1.98, 9.15] 
3.75 [2.05, 5.40] [.94, 7.76] 


13.28/20 
208.26/190 
193.30/190 


193 SNe 
binned flux 6 
unbinned flux 
unbinned mag. 


.35 
.36 
.36 


(Zmax = l-199) 

[.28, .41] [.23, .48] 
[.30, .42] [.25, .49] 
[.30, .42] [.25, .48] 


.660 
.656 
.662 


[.645, .674] [.631, 
[.645, .664] [.637, 
[.653, .670] [.645, 


.688] 
.672] 
.678] 


1.39 [1.08, 1.69] [.84, 2.10] 
1.44 [1.06, 1.82] (.80, 2.41] 
1.42 [1.13, 1.71] [.89, 2.08] 


4.95 [2.55, 7.34] [.88, 10.35] 
7.50 [4.39, 10.69] [2.57, 15.57] 
5.14 [2.88, 7.43] [1.41, 10.62] 


14.24/21 
209.42/191 
194.50/191 


194 SNe 




(yZmax — 1.755) 














binned flux 6 
unbinned flux 
unbinned mag. 


.40 

.38 
.38 


[.32, .48] [.25, .54] 
[.32, .43] [.27, .48] 
[.31, .44] [.26, .50] 


.661 
.654 
.662 


[.646, .675] [.632, 
[.644, .662] [.637, 
[.652, .669] [.645, 


.688] 
.670] 
.678] 


3.26 [1.76, 4.76] [1.02, 5.75] 
2.85 [1.88, 3.81][1.18, 4.69] 
2.48 [1.66, 3.54][1.22, 4.35] 


15.64 [6.22, 25.30] [.92, 30.53] 
14.78 [8.77, 20.58] [4.36, 25.83] 
10.60 [5.41, 17.81] [2.68, 21.80] 


14.5/22 
209.88/192 
194.70/192 


58 SNe 




(Zmax =0.863) 




Knop sample 

offsct-23.5 




px(-5 Z rn ax)/p(0) 


PX (Zmax) / p(0) 


Xmin/ N dof 


binned flux 6 
unbinned flux 
unbinned mag. 


.34 
.26 
.39 


[.25, .43] [.18, .52] 
[.18, .33] [.12, .41] 
[.30 .48] [.22, .56] 


.311 
.381 
.309 


[.210, .409] [.127, .509] 
[.331, .428] [.286, .474] 
[.256, .362] [.208, .412] 


.91 [.59, 1.22] [.38, 1.61] 
.85 [.67, 1.03] [.50, 1.21] 
1.18 [.89, 1.46] [.67, 1.82] 


5.92 [2.13, 9.73] [.58, 13.21] 

1.96 [.53, 3.55] [.09, 5.45] 
6.10 [2.47, 9.59] [.86, 12.57] 


7.2/12 
61.32/56 
53.90/56 



° Statistical error only, not including the contribution from the much larger SN la absolute magnitude error of a™ 1 ~ 0.05 (see Sec. 4). 
6 flux-averaged with Az = .05. c The number of degrees of freedom. 
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Fig. 1. — The marginalized probability distributions of O m , h, and the constant equation of state for the dark energy uiq. 



CMB+LSS+(192 To nry/Barris sample SNe, flux- averaging st at. [solid] & mag stat. [dot]; z m .,= 1.056) 
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12 3 
Px(0-5 z m J/p x (0) 



5 10 15 
Px(z m «)/Px(0) 



CMB+LSS+(194 To nry/Barris sample SNe, flux- averaging st at. [solid] & mag stat. [dot]; z m ..= l .755) 







2 4 6 8 
Px(0-5 z mM )/Px(0) 



CMB + LSS + (58 Knop sampl e SNe, flux -avera ging stat. [s olid] & mag stat. [dot]; z, 
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Fig. 2. — The marginalized probability distributions of Q m , h, and dimensionless dark energy density at z = z max /2 and z = z., 
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Fig. 3. — The dark energy density px(z) reconstructed from current SN la (Tonry/Barris sample and Knop sample), CMB (WMAP, CBI, 
ACBAR), and LSS (2dF) data. The densely (sparsely) shaded regions are the 68.3% (95%) C.L. of p x (z) for 192 Tonry/Barris sample SNe 
la (at z < 1.056). The heavy (light) lines indicate the 68.3% (95%) C.L. of the reconstructed px(z). The dot-dashed line indicates the 
cosmological constant model, px(z)/px(0) = 1. (a) Reconstructed px{z) using 192 (shaded regions), 193 (dotted lines), and 194 (dashed 
lines) SNe la from the Tonry/Barris sample (Tonry et al. 2003; Barris et al. 2003). (b) Reconstructed px{z) using 192 Tonry/Barris sample 
SNe la (shaded regions) and 58 Knop sample SNe la (dotted lines) (Knop et al. 2003). 



